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SUR 

L'ÉVALUATION APPROCHÉE 



DES AIRES PLANES 



Objet du présent mémoire. 

1 . Formules anciennes. Il existe un grand nombre de formules 
pour le calcul approché des aires planes. La plus ancienne, dont 
l'idée première est due à Côtes et que Newton a insérée dans 
les Principes (Livre III, lemme V), repose sur la substitution, 
à la courbe donnée, d'une parabole de degré n, ayant avec 
cette courbe (n -h 1) points communs. L'expression approchée 
de l'aire ne contient que les (n -h 1) ordonnées communes à la 
courbe et à la parabole et les distances qui les séparent. Gauss 
a prouvé (Werke,lll, 165-196, 202-206) qu'en choisissant 
convenablement ces distances, oii pouvait obtenir une formule 
rigoureusement exacte, dans le cas où la courbe donnée est une 
seconde parabole de degré (2n — 1) (*). Malheureusement, ni 
Gauss, ni aucun de ses continuateurs n'a cherché, que nous 
sachions, la valeur maxima de l'erreur que comporte sa for- 



(*) Voir dans la Nouvelle correspondance mathématique , t. VI, pp. 396-402, une 
démonstration très simple de ce théorème de Gauss, démonstration due à M. E. Catalan. 
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mule, dans le cas où la courbe dont on veut obtenir Taire appro- 
chée est une ligne continue quelconque, au lieu d'être une 
parabole de degré (2w — 1). 

La formule sommatoire de Maclaurin et d'Euler n'a pas cet 
inconvénient : on est parvenu à trouver la valeur de l'erreur 
maxima à laquelle son emploi peut conduire (*). Mais elle pré- 
sente d'autres désavantages. L'aire cherchée, pour une courbe 
y = fx,est exprimée, non seulement au moyen d'un certain 
nombre de valeurs de la fonction /", mais aussi des dérivées suc- 
cessives de cette fonction. En outre, l'erreur qu'elle comporte 
ne peut être calculée d'une manière approchée qu'en cherchant 
entre quelles limites varie une certaine dérivée de la même 
fonction. Par suite, la formule de Maclaurin et. d'Euler est, en 
pratique, d'un emploi extrêmement peu commode. 

On peut déduire de la formule de Maclaurin et d'Euler, direc- 
tement ou indirectement, un grand nombre de formules plus 
simples et plus pratiques : la formule des trapèzes, la (première) 
formule de Simpson, la formule de M. Catalan, la formule de 
Poncelet, la formule de M. Parmentier, etc. La (première) for- 
mule de Simpson, qui semble de beaucoup la plus exacte de 
toutes celles que nous venons de citer, peut s'obtenir aussi 
comme cas particulier de celles de Côtes, de Newton ou de 
Gauss. 

2. Objet de ce mémoire. La démonstration habituelle des for- 
mules de Poncelet et de M. Parmentier ne laisse rien à désirer 
au point de vue de la rigueur ou de la simplicité et elle conduit 
tout naturellement à une estimation approchée de Terreur com- 
mise en les employant. Mais il n'en est pas de même pour les 
autres formules, qui, au fond, sont établies d'une manière à peu 
près empirique. Récemment, en étudiant les écrits de M. le 
général Parmentier sur les quadratures approchées, nous nous 
sommes aperçu qu'on pouvait démontrer toutes les formules 



(*) Voir, par exemple, Houel, Calcul infinitésimal, U I, pp. 473-487. 
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énumérées plus haut avec la même facilité que celle de Poneelet, 
estimer l'erreur maxima correspondante et la représenter géo- 
métriquement, souvent de plusieurs manières. De plus, on peut 
trouver un grand nombre d'autres formules, dont une au moins 
(dueàM.Ch. Dupain), deux peut-être, semblent avoir une valeur 
pratique, à cause de leur simplicité et de leur exactitude. 

Nous nous proposons d'exposer ici, d'une manière systéma- 
tique, l'ensemble des recherches auxquelles nous avons été 
conduit. La partie la plus élémentaire de cette étude a déjà 
paru dans la Mathesis (t. I, pp. 17-22, 55-36, février et mars 
1881). 

3. Principaux résultats obtenus. Parmi les résultats obtenus , 
nous signalons spécialement, au point de vue pratique, l'estima- 
tion graphique de l'erreur que comportent les diverses formules 
et particulièrement celle des trapèzes. 

Si l'on représente par AD l'erreur maxima pour cette dernière 
formule, les erreurs analogues, pour les autres formules, sont 
respectivement : 

i 

- AD pour la formule de Poneelet ; 

- AD pour la première formule de Simpson, et pour 
5 

celle de M. Parmentier et la formule de 
M. Dupain (n° 7); 

-AD pour la seconde formule de Simpson et pour 
la formule de M. Catalan; 

5 

-AD pour les formules de Poneelet, de M. Parmen- 
tier et de Simpson, quand on les étend au 
cas d'un nombre pair d'ordonnées (n° 21); 

4 

-AD pour la formule nouvelle (n° 15). 

5 
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Ces limites de Terreur peuvent servir, dans un cas donné, à 
voir rapidement, sur la figure même, par une construction plus 
simple que celle que Poncelet a indiquée pour sa formule, quelle 
grandeur il convient de donner à h pour obtenir une approxima- 
tion déterminée d'avance. 

Une fois A choisi, les calculs qui servent à trouver une valeur 
approchée de S, par Tune quelconque de ces formules (celles de 
Poncelet et de M. Parmentier exceptées), donnent en même temps 
une nouvelle limite de l'erreur commise et cette nouvelle limite 
est inférieure à celle qui vient d'être indiquée, même pour la 
formule des trapèzes. 

Voici encore deux résultats sur lesquels nous attirons l'atten- 
tion du lecteur : 

1° Quand le nombre des ordonnées est pair, la meilleure for- 
mule pratique est celle de M. Catalan , qui est, pour ce cas, la 
transformée de celle de Simpson , et qui peut d'ailleurs être 
employée même si le nombre des ordonnées est impair (n° 21). 

2° La formule de M. Dupain (n° 7), presque aussi simple que 
celle de Poncelet, a une exactitude presque égale à celle de 
Simpson, et est aussi exacte que celle de M. Catalan. 

11 

Préliminaires. 

4L. Estimation de l'erreur commise dans les approximations. 
I. Si une quantité inconnue S est plus grande que m 9 plus petite 
que M, on peut dire que S = m, avec une erreur,, par défaut, 
plus petite que M — m, ou que S = M, avec une erreur, par 
excès, plus petite que M — m. 

II. Si l'on prend S = ^ (M -+- m), l'erreur commise est tout 
au plus égale à |(M — -m), mais on ignore quel en est le sens. 

III. Si l'on sait que S est plus rapproché de M que de m, 
S étant alors supérieur à | (M -t- m), on peut prendre pour S 
une valeur intermédiaire entre M et |- (M h- m), par exemple 
| (2M h- m). En posant S = ~ (M h- m), l'erreur sera par excès 
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si S est compris entre \ (M -h m) et § (2M -4- w?) , et elle sera 
moindre que 

1 1 1 

-(2M ■+- m) (M -h wî) = -(M —m). 

L'erreur sera par défaut si S est compris entre M et £ (2M -f-m), 
et elle sera moindre que 

4 1 

M (2M -*-m) = -(M — m). 

3 3 

Donc, en tout cas, si S est compris entre \ (M -h m) et M, et 
si Ton pose S=f (2M h- m), Terreur est moindre que \ (M — m). 

IV. En général, si p est une quantité comprise entre m et M, 
en faisant S = fx, on commet une erreur plus petite que la plus 
grande des deux quantités p — m, M — p.. 

5. Données et notations. Considérons (Fig. 1, 2 ou 5) une 
aire S, comprise entre un arc de courbe AL, une droite fixe al 
et les perpendiculaires Aa, Ll abaissées sur cette droite des 
extrémités de l'arc AL. Supposons, pour fixer les idées, que 
l'arc AL soit, dans toute son étendue, concave vers la droite 
al (*). Divisons la base al en un nombre pair de parties égales, 
en dix, par exemple. Par les points de subdivision, 6, c, d, e 9 f, 
g, h, i, k 9 élevons des perpendiculaires à al. Désignons ces per- 
pendiculaires ou ordonnées aA, 6B, .... /L, dont les extrémités 
A, B,...L sont sur la courbe, par y<\,y 2 > •••> Vm\ soit h la distance 
commune de deux quelconques de ces ordonnées. Nous repré- 
senterons par E, la somme des deux ordonnées extrêmes, par I, 
la somme des autres ordonnées de rang impair, par P, celle des 
ordonnées de rang pair. De cette manière, 

1 = ^4- y* + yi + y*> 

P s =yj + ï4 + y6+y8 + y^ 



{*) S'il était convexe, les quantités d, d', d",J t seraient négatives aiHieu d'être positives. 
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Nous posons ensuite 



y% + yio yi -+- y u 



2 2 

d"=%d~ d' } 
(J = 2P — E— 2Ï 



III 



Formules de Poncelet, de M. Parmentier et de M. Dupain. 

6. Première limite supérieure et première limite inférieure 
de S. I. Par les extrémités B, D, F, H, K des ordonnées de rang 
pair (Fig. 1), menons à la courbe des tangentes terminées enB 1 




Fig. 1. 

et B 2 , D 1 et \) 2 ,F { et F 2 , H { et H 2 , K, et K 2 , aux ordonnées 
voisines. La somme des trapèzes 

ttE^BaC, cDiD 2 e 5 eFjF^, </H|H 2 i, iK ( K,l, 
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est supérieure à S. Ces trapèzes ont pour mesure 

2% 2 , 2% 4 , 2% 6 , 2% 8 , 2% 10 . 
Donc, 

M = A (2y, -h 2y, 4- 2y 6 -^ 2y 8 -h 22/ JO ) = 2ftP. 

IL Le polygone aABDFHKL/ a une aire m inférieure à S. 
On a d'ailleurs 

m s==h U (#»+#) -KîM-y*) +(y*+y«)+(»6+»8)+(»8-+-yio)H--(yio-4-yii) 

= A [2y 8 -*- 2y 4 -h 2*/ 6 -4- 2 </ 8 ^- 2y 10 ] — A - (y 2 -f- y i0 ) — - (y i -h t/u) , 

ou, d'après les notations du numéro précédent, 

m = 2/*P — M 

III. Il est impossible de démontrer d'une manière générale 
que S est plus rapproché de M que de. m. Néanmoins, dans le cas 
de courbes à courbure peu prononcée et lorsque h est suffisam- 
ment petit, on reconnaît aisément, en faisant un tracé graphique, 
que l'on a souvent le droit, au point de vue pratique, de sup- 
poser S plus rapproché de M que de m. 

7. Formules de Poncelet, de M. Parmentier et de M. Dupain; 
interprétation géométrique de l'erreur. I. Puisque S est com- 
pris entre M et m, on peut, d'après le principe II du n° 4, 
prendre approximativement M = p, p étant la demi-somme de 
M et de m; l'erreur s sera inférieure à la demi-différence. On 
aura ainsi la formule de Poncelet : 

S =p (approximativement), 

4 i 

p = - (M + m) = 2AP — - hd, 

1 i 
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II. Si Ton prend p' = f (2M -+- m), en faisant 

S = p r (approximativement) , 

on aura la formule de M. Parmentier. D'ailleurs, 

1 

p'=2AP hd. 

3 

L'erreur maxima est égale à la plus grande des deux quantités 
M — p , ==~hd,et p' — m=\hd, donc, égale à f hd. 

Très souvent, d'après le principe III du n° 4, et la remarque 
faite au numéro précédent, l'erreur sera inférieure à { (M — m), 
c'est-à-dire à f hd. 

III. Nous verrons plus loin pourquoi M. Parmentier a choisi 
pour p' la valeur | (2M -4- m), plutôt que ~ (3M •+-' m) , ou toute 
autre valeur comprise entre M et { (M -+- m). Nous verrons en 
même temps qu'une formule préférable à celle de Poncelet et à 
celle de M. Parmentier est la suivante, due à M. Dupain, 

S = 7r 4 (approximativement), 

1 hd 1 , , \ hd 

7r { = p' = c 2hV hd— > 

1 2w-l 3 3 'In — 1 3 

2Ad i hd 

o In — \ 3 

2n étant le nombre des subdivisions h de la base al, donc, dix 
dans le cas de la figure. 

IV. La quantité d et, par suite , les produits hd, ^hd, etc., sont 
susceptibles d'une interprétation géométrique simple. Appelons 
P, Q, les points d'intersection des droites AL, BK avec l'ordonnée 
moyenne F/". On aura évidemment 

1 i 

Par conséquent, 
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Il est donc facile de tracer sur la figure même un petit rec- 
tangle PQQ'P' dont Faire est égale à hd. 

Par suite, on peut estimer, sur la figure même, avant d'entre- 
prendre aucun calcul, la grandeur des erreurs théoriques { M, 
| hd, etc. Si ces erreurs, pour une valeur de h choisie d'avance, 
étaient trop considérables pour le but que l'on se propose en 
cherchant S approximativement, on pourrait recourir à une 
valeur plus petite de h, et ne commencer à estimer numérique- 
ment p, p', ou n { , qu'avec la certitude de trouver ainsi une valeur 
suffisamment approchée de S (*). 



Formule des trapèzes et formule de Simpson. 
Première démonstration. 

8. Deuxième et troisième limite inférieure de S. I. Le polygone 
aACEGIL/ (Fig. 2) a une aire m' inférieure à S (**). On trouve, 
sans peine, 

m' = h [{y { ■+■ y z ) -h (y z 4- y s ) -+- (y, -h y 7 ) -*- (y 7 -t- y 9 ) -t- (y 9 -4- y {l )\ 

ou encore 

m' = h (E h- 21). 

IL Le polygone aABCDEFGHIKL/ (**■*). a aussi une aire t 
inférieure à S. L'aire t est donnée par la formule 

1 

ou 

■r==-A(E-f-2I +2P). 

(*) Ce paragraphe ne contient absolument rien de neuf et doit être regardé comme 
une introduction aux suivants. 

(**) On ne conçoit pas comment Poncelet, à qui appartient la considération des aires 
M et m, n'ait pas songé à ce polygone m' qui va nous donner si simplement la démonstra- 
tration de la formule de Simpson. 

(*-**) Non tracé sur la figure. 
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9. Formule des trapèzes; formule de Simpson. I. La figure 
ou les formules donnent immédiatement 

*=i(M + m'). 

Donc, d'après le principe H, n° 4, en posant (formule des 
trapèzes) 

S = t (approximativement), 

on commettra une erreur maxima plus petite que ~ (M — m'). 
Or, 

M — m' = A(2P — E — 21) = A& 

L'erreur maxima, par défaut, que comporte la formule des 
trapèzes, est donc plus petite que { hd. 






Fig. 2. 

IL Puisque S est compris entre M et t = | (M -h w'), on 
pourra prendre 

S = s (approximativement), 

5 étant égale à | (2M + m'). C'est la formule de Simpson. On a 
d'ailleurs 

S :=ifc(E-4-2I-*-4P), 
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et d'après 4, III, 



s < - (M — m) ou - M. 
5 5 

On peut encore écrire s et /$ sous la forme 

s= i (m -*- 2e), ■ M = 2 (M — ' t). 

3 

On verra plus bas pourquoi la valeur la plus convenable pour s, 
quantité qui doit être choisie entre M et m\ ou M et t, est celle 
qui est donnée ici. 

10. Interprétation géométrique de 8. Autres limites de l'erreur. 
I. On peut donner, de à aussi, une interprétation assez simple. 
Soient R, S, T les points où CI, DH, EG rencontrent ¥f. On 
aura : 

y»- + -y» :=R/ - y* + y« _ Sjr y^^-y^ ^_ Tf 

2 2 2 

Or, la valeur de 5 peut s'écrire successivement : 

<te (2y 8 + 2y, + 2t/ 6 -t- 2y 8 -h 2y 10 - y* - y M - 2y 5 - 2y s - 2*/ 7 - 2y 9 ) 

«sL 2 y» + y' t 2 yi-*-ys 2 yr«-y» 1 ç> y*+y*> g yi-*-y« 

' w 6 « o o o a 



= 2 



L — y» + y^ _ ( y*+yi y*+y* \ . + ( y*+ys y*+y> 
y* -+- yio \ _^ / y» -*- y«> yi ^y^l 



y3j^y«_ 

2 

Par conséquent, 

^ =s2 [ (F /L.T/)-(Tf-3/) + (8/--Rn-(RA- Q/) + (QA- P/)] 
= 2 [FT - TS + SR — RQ + QP], 

quantité que Ton peut construire sur la figure. On peut donc 
aussi représenter h$ par un petit rectangle. 

II. La construction graphique de $ n'est pas assez simple pour 
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permettre l'estimation rapide de l'erreur correspondant à une 
distance h des ordonnées, quand on emploie la formule des tra- 
pèzes ou celle de Simpson. Mais il est facile de prouver que à est 
inférieur à 2d et que, par suite, Terreur, dans la première, est 
moindre que hd, et que f hd 9 dans la seconde. 

Pour le prouver, remarquons d'abord que l'on peut écrire : 

j = 2 [PQ — (QR — RS) — (ST — TF)]. 

Or, on a PQ = d, et PQ > QR > RS > ST > TF. Pour établir 
ce dernier point, prouvons, par exemple, que PQ ou d est supé- 
rieur à d' ou QR. L'inégalité à établir, PQ > QR, devient succes- 
sivement : 

11 4 \ 

5 (2/2 + y*o) — - (y* + y a) > ^ (y 3 + v*) — ^ (y* + y™), 
\ \ 

2/2 + 2/io — -(2/1 -t-y u ) — -(y» + y 9 )>0, 
4 4 

y 2 — - (y* '+ y s ) + y w — ~ (y 9 + y«) > o. 



Soient, b { le point d'intersection de 6B avec AC, #*, celui de kli 
avec IL. On a 

4 4 

Donc l'inégalité à vérifier devient : 

B6 4 -+-Kfc 4 >0, 

et, sous cette forme, elle est évidente. 

Il est donc aussi facile d'estimer graphiquement une limite 
supérieure de l'erreur pour la formule des trapèzes ou pour celle 
de Simpson, que pour la formule de Poncelet ou pour celle de 
M. Parmentier. 
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Démonstration graphique simultanée des résultats précédents. 
Première extension au cas où le nombre des ordonnées 
est pair. 



11. Formule des trapèzes. On peut donner de la quantité hd 
une interprétation graphique différente de celle du n° 7, appli- 
cable même au cas où le nombre des ordonnées est pair. 




Fig. 3. 

Considérons, par exemple, une aire2=aABCDEF/ > (Fig. 3) 
ayant encore sa concavité dirigée vers la base al, et ses ordon- 
nées ak, 6B,...., Ff, distantes entre elles d'une même quantité h. 
Cette aire est comprisejentre le polygone inscrit T = aABCDEF/^ 
et le polygone circonscrit aB l B 2 D i D 2 E 1 E 2 /'. 

La différence A entre ces deux polygones est plus petite que 
la somme des triangles 

ABB d , BCB 2 , CDD l5 DED 2 , EFE 2 . 

Par le point B, menons BX, prolongement de CB, By, Bd, Bg, 
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Be 2 , B9 respectivement parallèles à DC, DjDDa, DE, EE 2 , EF. 

Les triangles énumérés tantôt sont égaux respectivement aux 

triangles 

ABB l5 BiBX, yB<?, 6Be, e 2 B f , 

indiqués sur la figure par des hachures. Leur somme est infé- 
rieure au triangle AB<p, qui est égal lui-même à la somme des 
deux triangles BPA, EQF, BP, EQ étant des parallèles à af. 
L'aire 2 est donc inférieure à Taire T augmentée de ces deux 
triangles, autrement dit, on a 

aABCDEF/'< aire 2 < aPBCDEQ/, 

résultat extrêmement simple et non encore signalé, que nous 
sachions (*). 

Si la courbe était continuellement croissante (Fig. 4) ou 
décroissante, de A à F, la différence A est encore inférieure à 
un triangle ABcp, dont le côté Bcp est parallèle à la dernière 
corde du polygone T. Ce triangle ABcp est égal à la différence 
des triangles BPA, EQF. Par suite T 4-ABP— EQF, est une 
limite supérieure de 2, qui est encore égale au polygone 
aPBCDEQ/i comme dans le cas précédent. 

Pour une figure donnée, il est évidemment plus facile de 
construire les deux triangles ABP, EQF que le petit rectangle 
PQQ'P' de Poncelet (Fig. 1). 

Dans le cas de la figure 3, on a d'ailleurs 

APQ = ^Hy 2 -2/0, DEF = - A (*,„_< -yj. 
Dans le cas de la figure 4, 



(*) Si les ordonnées n'étaient pas équidistantes, il est clair que l'erreur commise en 
prenant S=T est moindre qu'un triangle semblable à AB^? et ayant pour hauteur, 
au lieu de AP = /*, la distance des deux ordonnées consécutives de la courbe les plus 
éloignées l'une de l'autre. 
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Donc ABcp, somme des deux triangles dans le premier cas, 
différence dans le second, est égale à 

i 1 H i 1 

2 h (?/t — yO ■*- 2 * (y.-i — y») « * I g (y» ■+■ y-0 — 5 (y* -*- y»)]- 

quantité que Ton peut désigner par hDl 




Fig. 4. 

Appelons T Taire aPBCDEQ/, de sorte que T'=T + M). 

On aura : 

T = A (- y { h- 2/ 2 h- y, + . . - + */„_, -+- - </„) , 

T==zh (2 y " 4- ?/ 2 -h 2/ 3 4- • • • + 2/ n _ 4 -+- — 2/«— *)" 
On obtient donc une limite supérieure T' de Taire 2, en 

2 
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remplaçant dans Taire du polygone inscrit T, la première et la 
dernière ordonnée par la seconde et l'avant-dernière (*).. 

1£. Formules de Poncelet, de M. Parmentier et de Simpson. 
Ces formules n'existent que dans le cas où le nombre des ordon- 
nées est impair. Considérons donc seulement l'aire S—aABCDEe 
(Fig. 3), Le polygone m* inscrit de Poncelet est aABDEe. Pro- 
longeons DB jusqu'à sa rencontre en (3 avec &\. La différence 
M=M — m, entre le polygone circonscrit et le polygone inscrit 
est égale à la somme des quatre triangles 

ABB 1? BcA = B d Bj3, DjDjC,, = |3B£, DED 2 = âBe, 

ou au triangle ABs. L'erreur maxima, dans la formule de Pon- 
celet, est donc moindre que | ABs; dans celle de Parmentier, 
elle est inférieure à § ABs, et souvent à f ABs. 

On a vu, au n° 9, que dans la formule de Simpson, e<hd et 
A5=2 (M — t). Dans le cas actuel, d'après le n° 11, M—t est la 
somme des triangles compris à l'intérieur de ABs et marqués 
par des hachures. Donc A<î<2 ABs, |A<î<§ ABs et, par consé- 
quent, l'erreur s est moindre que § hd. 

18. Première extension au cas d'un nombre pair d'ordonnées.. 
Ce qui précède suggère un moyen d'étendre les formules de 
Simpson, de Poncelet et de M. Parmentier au cas où le nombre 
des ordonnées est pair. Soit à calculer approximativement l'aire 

2=S-t-<7, o- désignant l'aire curviligne eEFf. On a 

£<s<f-*-AB* 
trapèze eEFf = trapèze eEP/" < trapèze eEF/* -+- EE 2 F. 



(*) La formule des trapèzes (comme aussi celle de Simpson et les autres) est employée 
à l'estimation des volumes, des centres de gravité, des moments d'inertie, etc. Dans le 
jaugeage des navires, par exemple, y u ij 2 ,...., y n représentent les aires de sections équi- 
distantes et on peut les obtenir sur des coupes bien dessinées, au moyen du planimètre. 
On conçoit aisément qu'au moyen de deux lectures cet instrument donne même immé- 
diatement T. D'après la formule démontrée dans le texte, on peut de même trouver T'. 
Donc, au moyen de quatre lectures, on peut obtenir une limite supérieure et une limite 
inférieure de la capacité du navire. 
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Ajoutons, en remarquant que 

\ 

trapèze eEF/*= - h (?/„__, h- y n \ t -+- trapèze eEF/'= T, 

A 

ABs -+- EE 2 F ou ABf -*- e 2 B? < AB? ou hï). 
Il viendra : 

T < « + 2 & (y-* + ^ < T + ftD . ou T', 

Donc, en complétant la formule de Simpson, au moyen d'un 
trapèze inscrit, dans le cas d'un nombre pair d'ordonnées, l'er- 
reur est moindre que M); autrement dit, la valeur trouvée, 
comme 2, est comprise entre les polygones .aABCDEF/ 1 , 
ePBCDEQ/1 

Dans le cas des formules de Poncelet et de M. Parmentier, 
on prouvera de même que T ; surpasse p ou p f augmenté de 
1 h (2/n-i"t-2/«)> ma * s ces quantités peuvent être inférieures à T. 
Néanmoins, on établit encore aisément que les formules approxi- 
matives 

1 1 

2 — P -*- - h {y n . { h- y n ), 2 = p' + - {y n -i + y n ) , 

comportent une erreur moindre que M). On a, en effet, 

A 

S = p db une partie de - ABe 



i 

* = -h (y»., + t/ n ) -t- une partie de j? a B ? * 



Donc, en ajoutant, 



i 

S-H(7 = 2 = p -i — A (t/„_i-t-t/ n )± une partie de AB^. 



De même, 



4 

2 = p' -t- -h(y n _i -+- y n ) =h une partie de ABp. 

À 
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VI 



Nouvelle limite supérieure de S. Formule nouvelle. Troisième 

DÉMONSTRATION DE LA FORMULE DE SlMPSON. 

14. Nouvelle limite supérieure de S. I. Par les extrémités 
C,E,G, I des ordonnées de rang impair autres que les extrêmes, 
menons, à la courbe (Fig. .5), les tangentes terminées en C 1 et C 2 , 




Fig. 5. 

E,, et E 2 , G { et G 2 , I>i et I 2 , aux ordonnées voisines. Prolongeons 
en outre CB, jusqu'à sa rencontre, en X, avec aA prolongé et IK 
jusqu'à sa rencontre, en Y, avec IL prolongé. On aura : 

aX * cG = 26B, IY + H = 2&K, 



ou, successivement, 



aX +y 5 = 2y a , /Y + y. = 2y 10 , 
aX = 2y a — y 3 , /Y = 2y 40 — y 9 . 
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La somme M' des trapèzes 

aXB6, bCfiA d^Ej/; fGfiA «A*,- &KY/, 

est supérieure à S. Ces trapèzes ont pour mesure 

4 1 

-&(3y,— y 5 ), % 8 , 2%„ 2% 7 , 2% 9 , -h{Zy lQ — y 9 ). 

Donc, après quelques réductions, 

M'=fc(E-*-2I)-ï-A(y,+.yio— î/i — y M ) — * Lte»-*-»»-- gfoi+yio) ; 

ou encore, 

M' = h (E -h 21) -4- 2ftd — Ad' = h (E h- 21) -t- M". 

II. On peut observer que d <d" <2d. En effet, on a vu plus 
haut que d' est inférieur à d; donc d"=2d — d' est positif et 
supérieur à d. 

15. Formule nouvelle. I. On obtient une formule analogue à 
celle de Poncelet et où n'entrent que des ordonnées de rang 
impair, y 2 et y i0 exceptées, , en posant 

S = p" (approximativement), 

1 1 

p" = - (M' -h m') = A (E h- 21) -+- - Ad". 

On a, pour Terreur maxima: 

1 1 

e<-(M' — m') ou -Ad"<Ad. 

Jà 2à 

On trouve, de même, une formule analogue à celle de M. Par- 
mentier: 

S =p'" (approximativement), 

p"' = -(2M' -*- m') = A (E -*- 21) '+ -Ad", 
3 3 

2 4 

£ < - Ad" < - hd. 

3 «3 
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II. Nous verrons plus loin que les deux formules précédentes 
ne sont pas très exactes. Il n'en est pas de même de la suivante , 

S = tt 2 (approximativement), 

„ == p "> \ ? hd" = h (E + 21) -*- - hd" - hd", 

' -2»-*- 15. 3 2n-Kl3 

qui est nouvelle et semble avoir une valeur pratique et théorique 
analogue à la formule 8=714 , donnée plus haut. Malheureuse- 
ment la limite supérieure de Terreur est assez élevée, savoir : 

2 12 

- hd" — hd'\ 

5 2» +• 1 3 

%n étant le nombre des subdivisions de la base al. 

16. Troisième démonstration de la formule de Simpson. I. Des 
inégalités 

M > S, M' > S, 
on déduit 

2M -*- M' > 5S, - (2M -*- M') > S. 
Or, 

i(2M -*- M') = -A(E-*-2I + 4P) + -Ad" = s -i- -Ad". 
5 3 3 3 

Donc, 

1 

s 4- _ hd" > S. 

5 
II. De même, des inégalités 

m < S, £ < S, 
on tire 

1 

m -*- 2£ < 3S, - (m ■+- 2«) < S 5 

Or, 

\ 1 1 1 

- [m + 2*) = - A (E -+■ 21 -*- 4P) Ad = s Ad. 

5 5 5 3 
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Donc, 

i 

S> s hd. 

3 

III. Puisque S est compris entre les quantités 

1 \ 

s hd, s+-hd", 

3 3 

on peut poser 

S = s (approximativement), 

avec une erreur moindre que f hd" 9 qui est la plus grande des 
deux quantités ~ hd, f hd". 

IV. Nous avons donc la limite de Terreur que comporte la 
formule de Simpson sous trois formes différentes : 1° s<|^; 
2°s< \hd; 5°g<{Ad". 

La seconde est la limite la plus grande, mais la plus facile 
à construire. La première est au contraire la plus difficile à con- 
struire, mais elle est toujours inférieure à la seconde, et, proba- 
blement aussi, à la troisième, car 

hâ = M — m', hd" = W — m\ 

et, en faisant un tracé graphique, on reconnaît que, le plus sou- 
vent M est inférieur à M'. Enfin, la troisième, certainement plus 
petite que la seconde, est presque aussi facile à construire (*). 



VII 

Nouvelles démonstrations de la formule des trapèzes 
dans tous les cas. 

17. Cas d'un nombre impair d'ordonnées. I. Le procédé qui a 
servi au numéro précédent à trouver de nouvelles limites de S. et 



(*) On peut prouver que hd" est inférieur à Qhd, par un tracé graphique analogue à celui 
du § V, mais cela ne présente guère d'intérêt. 
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une seconde démonstration de la formule de Simpson, peut évi- 
demment donner d'autres limites et d'autres formules en grand 
nombre. Voici, par exemple, une nouvelle démonstration de la 
formule des trapèzes, c'est-à-dire une nouvelle estimation de 
Terreur que l'on commet en l'employant. 

II. Des inégalités 

M > S, M' > S, 
on déduit 

M-v-M'>2S, *(M-f-M')>S. 
Or, 

i(M H- M')=i h (E -h 21 -t- 2P) + l-hd" = t ^-hd". 

On a donc : 

\ 
t ~f- - hd" > S > t. 
2 

Par suite (4, 1), si l'on fait 

S == t (approximativement), 

on commet une erreur, par défaut, inférieure à { hd". 

III. Posons encore : 



1 

9 



t -+- (t -+- -hd" 



\ 

. t -*- 7 hd", 



S = t' (approximativement). 

D'après le principe II, l'erreur maxima que comporte cette 
formule des trapèzes corrigée est inférieure à { hd" et le sens 
en est inconnu. 

18. Cas d\m nombre pair d'ordonnées. Prolongeons la 
courbe AL (Fig. 6), jusqu'à une douzième ordonnée A^=j/ 12 , 
située à une distance h de U. Appelons <r l'aire curviligne /LaA, 
et cherchons une valeur approximative de S -ho*. Pour cela, 



2S - 



prolongeons, la corde KL jusqu'à sa rencontre, en £, avec ).A. 
L'aire ILÇl surpasse cf. On aura, d'après le n° 14, 

1 

tLfr=-h(3y ii — y l0 ). 

Ajoutant cette aire à M, qui est plus grande que S, il vient : 

h (Zy 2 + % 4 + 2?/ 6 -f- 2*/ 8 + -t/ 10 + - 2/iiJ > S ^- (T. 

On trouve, de même, 

(33 \ 

2 2/2 + 2 ^ "*" 2 ^ 5 + 22/7 "*" 22/9 "*" 2 H ^ S "*~ *' 

D'où, en ajoutant et divisant par 2 : 

4/7 3 3 7 \ 




Le premier membre de cette inégalité peut s'écrire : 

h (2» 1 -*- y« + y» ■+-'•••-*- y«> -*- y« -*- g y«) 
. ."*" \ 2 2 / 2 l 2 2 / 
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Posons : 
















h 


-»( 


1 


+ 2/2 


h- 2/ 3 


-H . . . 


-*- 2/io -*- 2/n 


+ 5^ 


d,= 


2 


2/n 


2/i 


2 


, d; 

2^- 


2/3 -*- 2/io 
2 
-d;. 


2/2 ■+• 2/n 
2 



On pourra écrire là dernière inégalité : 

i 

h -4- - Adï >S + <r. 

La quantité ^ est la mesure de la somme des trapèzes intérieurs 
inscrits dans S-4-ff. On a donc : 

h < S H- (T. 

On déduit des inégalités précédentes, ces deux conséquences : 
1° Si l'on fait 

S h- a = t { (approximativement), 

on commet une erreur, par défaut, inférieure à | Ad/. 2° Si Ton 

fait 

1 

S -t- o- = < 4 -i- - Ad/, 
4 

on commet une erreur de sens inconnu, inférieure à | Ad'/. 

f». Démonstration analytique de la formule des trapèzes 
dans tous les cas, sous la forme donnée au § V. I. Considérons 
l'aire- 2 d'une courbe dont la concavité est toujours tournée dans 
le même sens et dont on connaît n ordonnées équidistantes y { , 
y 2 , 2/ 3v .. ? y n _ i9 y n9 leur distance commune étant A. Posons: 



"»G 



2 2/i ■«- 2/2 h- 2/3 ' -*- • • * -»- 2/»-2 -h. 2/n-i -*- ^2/n 

1 1 4 1 

D = 5 (2/2 -*- 2/»-i) — 5 (yi + 2/n), D' = - (y 5 + y B _,) — - (y, -h y n _,). 
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Les résultats des deux numéros précédents peuvent s'écrire : 

T<2 <T-h AD — -hW; 

ce qui prouve que la formule des trapèzes 2= T (approximati- 
vement) comporte une erreur par défaut moindre que AD — | hD'. 
II. On remarquera que D surpasse D', comme il est facile de 
le prouver (n° 10). On a donc v à fortiori, 

T < I < T -+- /*D, 

ce qui équivaut au résultat démontré graphiquement au § V. 
En voici une démonstration analytique qui ne diffère pas essen- 
tiellement de celle du n° .11 (Fig. 1, à compléter par le lecteur). 
Ajoutons à un premier trapèze analogue à aBjB^c, compris 
entre y { et y z , des trapèzes analogues à ILÇk, compris entre ?/ 3 
et 2/4-» y il et 2/ s , etc. La somme de ces trapèzes sera : 

h ^ 2 ^(^ 3 _i^H.g^_I^H- ... +(^-5 y^)] 
= h \\y* -*- y's •+- y* h — +.y«-3 -*- y«-> -*- ^ y«-i = t 

ce qu'il fallait démontrer. 



VIII 

Nouvelle extension des formules de Simpson, de Poncelet, etc., 
au cas d'un nombre pair d'ordonnées. Formule de M. Catalan. 

20.' Estimation approchée de l'aire g. I. L'aire a étant plus 
grande que le trapèze ILaI (Fig. 6), on a 

1 

* > ^ h to« + »«)■ 



AD: 
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D'après le n° 18, a* est inférieur à /L|>, ou 

1 

- h (3y M — 2/io) > ff. 

Toutes les valeurs intermédiaires entre les deux limites ici indi- 
quées peuvent servir à estimer approximativement or. Pour des 
raisons que nous indiquerons plus bas, l'a valeur la plus conve- 
nable est 

1 

o- = - (5/LAX ■+- /L?A) (approximativement), 
6 

c'est-à-dire, toute réduction faite, 

1 

M 

L'erreur maxima est inférieure à 

5 5 

g (/LÇx — ;LA>) = — h {— y l0 ■+- 2y„ — y„), 

d'après le principe IV du n° 4. 

II. On peut trouver aisément une aire v inférieure au trapèze 
IaLI et, par suite, à a et telle que 

^(2/LSx -h v)==~(5/La> h- /L?/). 
On trouve sans peine : 

On obtient géométriquement cette limite inférieure de <r, limite 
au reste toute artificielle, en menant a/ 2 parallèle à KL, L 2 \ 
parallèle à Ka. Alors v==ll 2 l 2 l. 

£1. Extension des formules de Simpson, de Poncelet, etc., au 
cas d'un nombre pair d'ordonnées. Formule de M. Catalan, 



comme extension de celle de Simpson. L Au lieu d'ajouter aux 
quantités s, p, p' (et aussi à n i9 p",p'", îr 2 ), Taire du trapèze 
ILa\ pour avoir une valeur approchée de S -ho-, comme on Ta 
indiqué au n° 13, il vaut mieux y ajouter la valeur approximative, 

1 

— h {—y io -*- Sy H + Sy«), 

trouvée ci-dessus. L'erreur commise dans l'estimation de S-+-c 
est égale à la somme des erreurs commises en cherchant séparé- 
ment S et c\ On prouve, sans peine, comme au n° 13, que cette 
erreur, dans le cas actuel, est inférieure, pour toutes les formules 
en question, à f AD. 

IL La formule de Simpson devient ainsi : 

1 

S "*" 12 k ^~~ yi0 "*" Sy " "*" ^ 

\ l 15 „ 5 

= - h Uji -+- 4y a -H- 2y 5 h h — #«> -h ay« ■*-£#« 

On trouve de même, pour une valeur approchée de la même 
aire : 

1/5 15 \ 

En prenant la moyenne de ces deux valeurs, on obtient l'expres- 
sion symétrique : 

„ f /3 7 23 ■ 23 7 3 

in i i i ri il 

*ï h h(y^y^-^y^y^\--'^ 11 ^^y^—^^y^j 

i i 

— L -*- - hd { hd[. 

4 12 
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La formule approximative, 

S = C, 

est due à M. Catalan. Nous allons la démontrer par une méthode 
au fond identique à la précédente, mais qui ne suppose pas que 
le nombre des ordonnées soit pair. 



Formule de M. Catalan. 

22. Formule de M. Catalan. I. Nous avons trouvé, aux n° 8 \\ 
et 19, que Taire 

, /3 3 

M 2 = h 1 5 2/2 ■+■ 2/3 ■+- • • • -*- 2/«-3 + y«-a -*- -y n -i 

dans le cas d'une courbe tournant en concavité vers Taxe al, 
supérieure à Taire 1 de cette courbe, dont on donne les n 
ordonnées équidistantes i/ 1? y 2 v> î/n- 
On peut regarder M 2 comme la somme des trapèzes ayant 

pour aire 

1 1 

2% 2 , - h (3y 3 — y*) , • • • , - h (5y n _ { — y n _ a ), . . (1 ) 

dont le premier seul est un trapèze circonscrit de Poncelet, et 
dont les autres sont analogues à la limite supérieure de o* consi- 
dérée au n° 20. On peut aussi regarder M 2 , comme la somme 
des trapèzes ayant pour aire 

1 1 1 

- h (5y a — y z \ -A(3y 8 — y t ), ..., -A(3y^ a — y,,^), 2%,^. (2) 

II. D'après le n° 20, la somme des trapèzes qui ont pour 
mesure 



(3: 



■a (yi-*- ys), ^ (^ y* — ïs + ^y* 
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dont le premier est un trapèze inscrit à 2, les autres sont ana- 
logues à la limite inférieure v de <r considérée au n° 20, est infé- 
rieure à 2. Cette somme est 

/ 3 3 i 5 \ 

* \ yi + 4^ 2 + 4 ^ 3 "*" ^ 4 + ^ 5 + " ' "*" y - f "*" 4 y "- 4 + 4 y "j" 

Pour la même raison, les aires 






(*) 



ont une somme, 

(51 3 3 \ 

inférieure à 2. 

Il en sera de même de la moyenne de ces deux sommes: 

, /9 \ 7 7 4 9 \ 

III. Nous avons dit, à propos de la formule de Simpson et au 
n° 20, et nous prouverons ultérieurement, qu'en général, la 
valeur la plus convenable à choisir pour chaque trapèze curvi- 
ligne composant 2 est égale aux f du trapèze appartenant à la 
suite (1) ou (2), augmenté de f du trapèze correspondant de la 
série (3) ou (4). Donc, si Ton pose 

i 

C=-(2M 2 -*- m 9 ), 

la valeur la plus convenable pour 2 sera donnée par la formule 
de M. Catalan, 

S = C (approximativement). 



— 52 
On trouve d'ailleurs : 



M 2 =T -h AD, m 8 = T — - AD — . - AD' , 

4 4 

G = T -h - AD — AD'. 



IV. Les principes généraux du n° 4 conduisent à une esti- 
mation de l'erreur beaucoup trop grande. On trouve s<C— m 2 , 

3 4 

G — w 2 = -AD ■+ -AD': 
2 6 

mais il est facile de trouver une limite moins élevée. Puisque 
U>D', on a C et I, compris d'après le n° 19, entre les quantités 

T, T + AD — -AD', 

Donc l'erreur commise, en posant 2=C, est moindre que la plus 
grande des quantités 

C — T = - AD — — AD', T + AD — i AD f — C =- AD — — AD', 
4 12 ' 2 - 4 12 

c'est-à-dire évidemment la seconde. En pratique, on peut dire 
que s est moindre que f AD. 



Seconde formule de Simpson. 

23. Élément de la seconde formule de Simpson. I. Considé- 
rons la partie aècdDCBA de l'aire S, et enfermons-la entre deux 
limites. 

Cette aire est plus petite que la somme des trapèzes 

«6BX, bCfiA 
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dont il est question au n° 14. L'aire de ceux-ci est 



ft 8 *-*-KM yi -./fcggLA 



On peut remarquer que le second de ces trapèzes a la même 
aire que le trapèze 6BD'd, D' étant le point d'intersection de BC 
avec Dd. On a donc : 

aire aXU'd=h ^^^\ 

2 

résultat qu'il est facile d'établir directement. 




Fig. 7. 

II. L'aire courbe abcdDCBA est supérieure à l'aire polygo- 
nale oôcdDCBA, qui a pour mesure 

1 

2 h (ï* "*" 2 »« + 2 2/3 h- »*). 

III. Nous verrons ultérieurement que la valeur la plus con- 
venable pour l'aire courbe aABCDd est 

\ n 5 -] 5 
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L'erreur commise est plus petite que la plus grande des deux 
différences, 

1 3 

- h (3y, -+- 3y B ) — g A [y* -*- 5y a + y 8 •+- y 4 ] 

3 1 

- h (y t + 3^ 2 + 3y, + y l ) — -h {y, h- 8y, +■ 2*/ 3 -+- y«) 



= fc 



i \ ■ \- \ 1 



c'est-à-dire que la première. Cette quantité est facile à inter- 
préter géométriquement. 

24. Seconde formule de Simpson. I. Supposons de même que 
Ton ait à chercher Taire 2 d'une courbe dont la concavité est 
tournée vers Paxe des x et dont on connaît (3n -h 1) ordonnées 
équidistantes 2/j, y 2 , ?/ 3 ,..., j/,^. On a > comme limite supérieure, 

3 

- h (% 2 + 2y 3 +■ 2y b h- 2y 6 -h 2y g + 2y 9 -*- ». + 2y5»-i ■+■ 2y 8B ) , 

4 

et, comme limite inférieure, 

-h{y { + fyz + 2y 8 -h . . . -h 2y 8fI -h y^). 

II. On prendra pour valeur approximative de Paire, le quart 
de la première expression plus les trois quarts de la seconde, 
c'est-à-dire : 

3 

o 

La formule 

S==s' (approximativement) 

est la seconde formule de Simpson. L'erreur maxima commise 
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est plus petite que la différence de s' avec la limite supérieure. 
On a donc 

3 P -, 

o 

III. On peut obtenir une autre limite de Terreur, plus facile 
à calculer à priori, par le procédé qui a servi plus haut pour 
la formule des trapèzes, au § V. On trouve immédiatement 
s <f M), D étant donné par la formule 

1 1 

D = g (y* + y ^) — ^ (2/l + y*.+i) > 

ou représenté par la somme ou la différence de deux triangles 
analogues à ABP, EQF de la figure 3. 



Formule de Weddle. 

25. Limite supérieure et limite inférieure. Considérons une 
courbe ABCDEFG concave vers Taxe ag, et sept ordonnées ka, 
B6,. e .., Gg, équidistantes d'une quantité h, que nous désignons 
encore par j/i,^,..., y v Estimons l'aire S 1 =aABCDEFG#. 

Par les extrémités B, D, F des ordonnées de rang pair, me- 
nons encore les tangentes B (1 B 2 , D d D 2 , F/jFg limitées aux ordon- 
nées voisines. La somme M { des trapèzes 

aB^c, cDjy^e, éF^g, 

est supérieure à Sj. Or, 

M 4 =2% a -i-y 4 -t-y 6 ). 

Menons les cordes AC, CE, EG. Le polygone aACEG# est 
inférieur à S^ et a pour mesure 
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Menons aussi les cordes AD, DG. Le trapèze aXDGg esl aussi 
inférieur à S/j et a pour mesure 

Des inégalités 

on déduit 

3mj -+• 2m 2 < oS 4 , - (3?^ -f- 2m 2 ) < S 4 . 
5 

Or, 

i 3 

m 3= - (3^ -t- 2wi a ) = - h (2î/ 4 -t- 2î/ 3 h- 2j/4 -*- 2.V& +• %7)« 

26. Formule de Weddle. L aire S 1 étant comprise entre M { 
et w 3 , on peut prendre pour valeur approximative de cette aire 




toute valeur comprise entre M 1 et m 3 . La valeur la plus conve- 
nable est, comme on le verra plus loin, 

\ 3 

io= - (3M^ m,) = — h(y { •*- Sy 8 + y, +. 6y< -+- y, h- 5t/ 8 -*- y,). 
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Nous poserons donc 

S, = w (approximativement), 

ce qui est la formule de Weddle. 

L'erreur commise est plus petite que 

T (3M 4 -t- m 3 ) — wt 5 = 7 (Mi — m 8 ) = 5(M 1 — w). 

4 4 

On peut déduire de ce qui précède une formule servant à 
trouver les aires dont la base est divisée en 6n parties égales. 



XII 



Exemples numériques (*). 

%7. Premier exemple. Soit à chercher Taire de la courbe dont 
l'ordonnée est égale au logarithme népérien de x, diminué de 3, 
x variant de 1000 à 1006 et h étant égal à l'unité. On trouve, en 
dix-billionièmes (l'aire exacte étant donnée par le calcul intégral) 



2/i = 0, 

t/ 2 = 43 407 74.79, 
y 5 === 86 77215.31, 
y, = 130 093 30.20, 
2/5=173 37128.09, 
y, = 216 606 17.57, 
y 7 = 259 79807.20. 

d =10793, 
ef= 6476, 
d"= 18110, 
6 =12951. 



S =78017128. 

s = S, . 

s' = S, . 

w = S, . 

c =s, . 

t = S — 2158 . 
p =S— 1080, 
p' =S + 720, . 

7T 4 = S. ... 

p" = S — 1079, 
p'"=S 4- 1439, 

> 2 == S . . . . 



« <4317 
e < 6475 
f < 12951 
fi < 4396 
fi < 6476 
£ < 5396 
s < 7195 
e < 6476 
fi < 7555 
e < 10073 
e < 8634 



(*) Les exemples traités ici ont exigé de nombreux calculs, pour le premier à douze 
décimales, pour les autres à huit décimales et avec deux ou trois groupes de données dis- 
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Les calculs ont été faits en allant jusqu'aux trillionièmes, c'est-à- 
dire, en prenant deux décimales de plus que celles qui sont con- 
servées. La première formule de Simpson et celle de Weddle 
donnent un résultat exact jusqu'aux trillionièmes. Même si l'on 
n'emploie que trois ordonnées, pour la première formule de 
Simpson, c'est-à-dire si l'on pose s=y { h- 4î/ 4 -4- y 7 , la différence 
S — s est moindre qu'un cent-billionième. 

Quand les erreurs théoriques peuvent être- calculées de plu- 
sieurs manières, nous n'avons inscrit dans la dernière colonne 
que l'erreur maxima la plus faible. 

28. Autres exemples. Dans le cas précédent, p, p', ^ sont 
calculés au moyen de cinq ordonnées seulement (?/ 3 et y% étant 
exclus), p", p"\ 7r 2 au moyen de 6 (y A étant exclu); les autres 
valeurs ont été trouvées au moyen des sept ordonnées. Voici 
trois autres exemples, où, pour chaque formule, on a employé 
seulement sept ordonnées : /& = £ pour les S premières, h===~ 
pour les trois suivantes ; enfin h=± pour les dernières. Tous les 
résultats sont donnés en millionièmes. 



Courbe I. 



Courbe II. 



Courbe III. 



Aire exacte 
Simpson I . 
Simpson II. 
Weddle . 
Catalan. . 
Trapèzes . 
Poncelet . 
Parmentier 
Dupain. . 



. S =693147 S =785398 S =4913223 

.s =S-4-23 s =S s = S — 3 

.s' =S-h48 s' =S-2 s' = S — 1 

. w = S-*-2 w =S-1 w =S 

. C = S-v-66 C =S-45 C =S — 10 

. t =S-h1730 t =S-1157 t =S — 134 

,p =S-f-376 p =S-231 p = S— 327 

. p' = S — 163 p' =S-M24 p' =S 4-102 

. ttj = S — 43 vi =S-*-45 tt 4 =S-t- 7 

p".= S4-603 p" =S-361 p"=S— 468 

p"'=S — 488 p'"=S-i-387 p'"=S-*-377 



Formule nouvelle . jt 2 = S — 3 



7T 2 =S-*-55 7T 2 =S-h2 



tinctes. Il se pourrait donc qu'il se fût glissé quelques erreurs dans nos résultats, mais 
elles ne sont pas bien considérables, pensons-nous. La courbe I a sa convexité tournée vers 
l'axe al, autrement dit, cl, cl', d", $ sont négatifs. 
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La courbe I est une hyperbole équilatère, de sorte que l'or- 
donnée est égale à [1 : oc], et x varie de 1 à 2. Dans la courbe II, 
l'ordonnée est égale à [1 : (1 + x 2 )], x variant depuis jusque 1. 
La courbe III est une circonférence de rayon égal à 2; on es- 
time l'aire comprise entre cette circonférence, un premier rayon, 
un second rayon perpendiculaire au premier et une parallèle au 
premier par le milieu du second. On peut prouver, par la géo- 
métrie élémentaire, que les aires de ces trois courbes sont pour I, 
le logarithme népérien de 2, pour II, \^ y pour III, |tt + 1 1/5. 

29. Exemples dans le cas d'un nombre pair d'ordonnées. 
Voici, pour une partie des formules précédentes, les résultats 
que l'on trouve pour les courbes I et II, en se servant de huit 
ordonnées. Pour les formules de Simpson et de Weddle, on les 
applique à l'aire comprise entre la première et la septième or- 
donnée; on ajoute ensuite au résultat trouvé une valeur de <r 
calculée comme il est dit au n° 20. 



Aire exacte 
Simpson I . 
Simpson II. 
Weddle. . 
Catalan . . 
Trapèzes . 



Courbe I. 

S =693147 

s = S-*- 4 
s' = S -h 18 
w==S-~ 8 
C = S + 59 
T=S + 1272 



Courbe II. 

S =785398 
s =S-*- 4 
s'=S-t- 16 
w = S-h 13 
C = S — 22 
t =S— 850 



XIII 



Choix d'une valeur (i intermédiaire entre une limite supérieure 

ET UNE LIMITE INFÉRIEURE DE S, COMME VALEUR APPROCHÉE DE S. 

30. Principe. Dans les paragraphes précédents, nous avons 
enfermé l'axe S à mesurer, entre celles de divers polygones 
M, M', etc., plus grands que S, et d'autres polygones plus petits 
m, m\ etc. Appelons, en général, L la limite supérieure de S, /la 



- 40 - 

limite inférieure. Ces polygones L, / sont des sommes de tra- 
pèzes ayant pour hauteur commune une longueur h, aussi petite 
qu'on le veut. Les différences L — S, S — /, positives pour les 
courbes à concavité tournée vers la base de Taire à mesurer, 
négatives dans le cas contraire, ont l'une et l'autre pour limite 0, 
quand h décroît indéfiniment. 

Soit R, le rapport inconnu, toujours positif des différences 
L — S, S — /. De la relation 

L — S 

7; r^R. 



on déduit 

L -+■ RI 

S = . 

1 + R 

Lorsque h décroît, en général, le rapport R varie. Dans beau- 
coup de cas, on peut en déterminer la limite r, pour /i=0, sans 
connaître la valeur exacte de S. Supposons que nous prenions, 
pour valeur approchée de S, la quantité p, définie par la relation 

L + rl 



i -4-r 

L'erreur réelle e ===== f/. — S sera donnée par la formule 

_L + rl L -+- R/ __ (L — /) (R — r) 
t_ "l -*-r ~~~ 1 -+- R "~(1 -4-r)(l -+-R)' 

Le dénominateur de s est supérieur à l'unité; le numérateur 
est le produit de deux facteurs dont chacun tend vers zéro en 
même temps que h, puisque lim L = lim/=S, 1imR=»r. 
Si, dans l'expression de p, on mettait au lieu de r une autre quan- 
tité positive, Terreur aurait pour numérateur le produit de 
deux quantités dont une seule, L — /, aurait pour limite zéro. 

Donc, dans les cas dont il s'agit, pour h suffisamment petit, 
il est préférable de prendre pour S la valeur f/. donnée plus haut. 
Dans les autres cas d'ailleurs, en général, on n'a aucune raison 
pour ne pas faire de même. 
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L'erreur s, comme on l'a dit au n° 4, IV, est plus petite que 
la plus grande des deux différences 




Si r > 1, s .<. L — p; si r < 1, s < ■ /x — /. 

SI. Cas ow le principe précèdent est applicable. Il existe un 
cas assez général, comme on le verra au § suivant, où Ton peut 
trouver r aisément et voir que p. est la valeur approchée de S la 
plus convenable. C'est celui où la courbe considérée est continue, 
où l'inclinaison de ses tangentes et sa courbure varient d'une 
manière continue et où, enfin, la courbure ne varie guère, soit 
pour la courbe entière, soit pour les diverses, parties dans les- 
quelles on est amené à la décomposer. 

Ànalytiquement, si ?/.= F'x est l'équation de la courbe, 
F'x, F'x, ¥ f "x doivent être des fonctions continues et F'"x varier 
peu entre les valeurs extrêmes de x que l'on considère. 

On observera que souvent l'on pourra vérifier rapidement ces 
conditions, soit au moyen d'un tracé graphique, soit au moyen 
de calculs assez simples, quand F'x ne sera pas une fonction 
trop compliquée. 

Dans le cas où l'on ne parvient pas à voir si ces conditions 
sont remplies, et où aucune propriété spéciale de la courbe 
ne force à choisir pour \x une valeur autre que celle qui est 
donnée plus haut, le mieux évidemment est de s'en tenir à cette 
valeur, qui a, sur toute autre, l'avantage de se justifier dans le cas 
général indiqué ci-dessus. 
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XIT 



CALCUL DE r. 



32. Cas cte trapèzes simples. I. Considérons d'abord le tra- 
pèze curviligne aABb (Fig. 9). Appelons la première ordonnée 
F'(#), la dernière F '(oc -h A). Ce trapèze a pour mesure 

F (a -h A) — F (a) = AF' -t- i A 2 F"-t- i AT'"^,), 

F', F" étant écrits, pour abréger, à la place de F'(ac), F"(ac), et oc? 1 
étant une valeur intermédiaire entre x et oc -h A. 




Fig. 9. 

Par le point A menons la tangente AA 2 terminée en A 2 à B6 
prolongé. On aura 

6A 2 = aA -h AF" (x) = F' h- AF". 

Par suite, le trapèze «AA 2 6, circonscrit au trapèze curviligne 
aABb, a pour mesure 



i h[F + (F' + AF")] = AF' h- * A'F". 
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Le trapèze inscrit aABô a pour mesure 

2 4 



^ A[F'(a;) •+- F'(x -*- A)] = AF' -h - A 2 F"-+- - A 3 F'"(x 2 ), 



a; 2 étant compris entre xetx-h-h. 
Dans le cas actuel, 

L — S 6 v ' 

R = 



s * I F "'( Xl )_i F '''(*o 

Par suite, si h tend vers zéro, et, par suite x { et x 2 vers x, on 
a lim R = r = 2. De plus, si F"'(x) varie peu de x à x -h A, 
R a aussi une valeur voisine de 2. 

On trouve le même résultat si l'on remplace le trapèze cir- 
conscrit aAA 2 6 par le trapèze aB^Bb, B d B étant la tangente 
en B, prolongée jusqu'à sa rencontre avec l'ordonnée ak. 

II. Soit 6 1 le point d'intersection de la corde AC avec 6B. 
On aura 

66 1 = - (aA -*- cC) = z [F' (x) -h F' (x -+ 2fe)] 

== -[2F' -h 2AF" -h 2A 2 F'"(x 3 )l = F' h- AF" -h A 2 F'"(x 3 ), 

x 3 étant compris entre x et x -h 2A. L'aire du trapèze akb { b 
sera 

-/i[F'+(F , +AF"+A 2 F'"(^))]=/iF'-+- - A 2 F" + i A 3 F'"(x 3 ). 

Si l'on prend cette aire pour limite inférieure, dans la valeur 
de R, il vient : 

L-8 -\ V "'M 

Q=5 — 



6 . v ' 2 ; 



valeur peu différente de ^, si F'"(œ) varie peu de x à x 4- 2A. 
Si A tend vers 0, on trouve lim R = r = f . 

On arrive au même résultat encore, si l'on remplace le tra- 
pèze circonscrit aAA 2 6 par aB d B6. 



III. Soit X le point d'intersection de CB avec Xa. On a, 
comme on l'a vu au n° 14, 

aX = 26B-- cC = 2F'(# -h A) — F'(x -t- 2A). 

Le trapèze aXBôa pour mesure 

-A[5F'(a;-4- A) — F'(ai+2*)] 

= I A pp h- 3/>F" -4- - A»F'"(ic«) — F' — 2/jF" — 2tfF'" (x s )] 

= AF' -*- - A 2 F" + j-A'Ff"(a? 4 )— A 5 F"(a?,), 

x 4 étant compris entre x et x h- A, % entre x et ac -+- 2A. 

Si Ton prend cette aire pour limite supérieure, dans la valeur 
de R, et akBb pour limite inférieure, il vient 

L-S ^'"K)-F'"W-^'"W 
R= 



1 If-^-If"^^ 



valeur voisine de 5, si F" '(oc) varie peu de x à x -+• 2A. Si A tend 
vers 0, on trouve lim R = r =5. 

33. Somme de trapèzes. I, Dans l'élément de la seconde for- 
mule de Simpson, faisons y l = F' (se), ?/ 2 = F' (x -h A), etc. 
L aire de cet élément sera 

F (a + 5A) -¥x = 3AF h- ~ A 2 F"-4- ? A 5 F'" (a? f ), 

x 6 étant compris entre x et x -4- 3A. 
La limite supérieure est 

-A[F'(a + A)n-F'(x4-2A)] 

= 5 A [2F' -h 3AF" -4- A 2 F'"(a; 7 ) + 2A 2 F'" (a;,)] 

ji 

=3ftF' + - A 2 F" -+- - fe 3 F'"(a; 7 ) + 3ft 3 F'"(x 8 ), 
2 4 

% étant compris entre x et x -4- A, x 8 entre x et x •+- 2A. 
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La limite inférieure est 

- h[F'x h- W{x +• h) 4- 2F' (x -h 2ft) -+• F'(x ■+- 3 A)] 

= 3AF' -4- -ft s F"-+- ift 8 F"(x 9 ) -+ 2ft 8 F"'(a? 10 )'-*- - A 3 F'"(xii), 
2 2 4 

x 9 étant compris entre x et œ ■+ A, x 10 entre a? et x -+- 2/*, sc^ 
entre x et x -*- 3h. 
Par suite, 

jF'"(a? 7 ) + 3F'"(* 8 )-^F'"(ar 6 ) 

R = ■ : ; : , 

~ F"> 6 ) - i F'" (x 9 ) - 2F'" (* 10 ) - - F'"(x ld ) 

valeur qui a pour limite r = 3, et qui est d'ailleurs voisine de 
3, si F'" x varie peu de oc à x -f- 3A. 

II. Des calculs absolument analogues prouvent que, dans le 
cas de l'élément Weddle, r = | et que R a une valeur voisine 
de |, pourvu que F f " varie peu de x à x -+- 6h. 



Application aux diverses formules précédentes. 

34. Formules de Simpson, de Weddle, et de M. Catalan. La 
première valeur de r trouvée au n° 32 prouve que la formule de 
Simpson, sous la forme 

S = s (approximativement), 

M -*- 2J i 

s== = - M -4-2*)» 

M ayant le sens indiqué au n° 6, est bien établie d'après le prin- 
cipe du § XIII. On voit de plus qu'elle sera probablement très 
exacte, si F'" (x) varie peu dans chaque intervalle h. 
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La seconde valeur de r (n° 32) conduit aussi à la valeur s, 

sous la forme 

i 
M -f- -m' 

2 \ 

s = = - (2M -h m'), 

4 -f-- 
2 

m' ayant le même sens qu'au n° 8. On ne voit pas aussi bien, au 
moyen de cette valeur de r, pourquoi la formule de Simpson 
est généralement si exacte; car, pour établir que R est voisin de 
r = {, on doit supposer ici que F'" (x) varie peu dans chaque 
intervalle 2A. 

La formule des trapèzes n'est pas aussi rationnelle que la for- 
mule de Simpson, puisqu'elle suppose que R, dans ce second cas, 
est égal à 1, valeur qui n'est pas voisine de f . 

La troisième valeur de r, combinée avec les considérations du 
n° 20, II, conduit naturellement à la formule de M. Catalan; la 
quatrième à la seconde formule de Simpson, la cinquième à 
celle de Weddle. On admet que F"' varie peu, dans un inter- 
valle %h pour la première, 3h pour la deuxième, &h pour la 
troisième. 

35. Formule de Poncelet, de Parmentier , de Dupain, etc. 
I. Pour une aire divisée en 2w trapèzes curvilignes, les limites 
M et m données au n° 6 sont telles que dans le premier et le 
dernier de ces trapèzes, la valeur de R est voisine de 2, pourvu 
que -F'" (x) varie peu dans un intervalle h; la valeur de R est 
voisine de |, pour les (2n — 2) autres trapèzes curvilignes, 
pourvu que F"' (x) varie peu dans un intervalle 2/i. 

On reconnaît aisément qu'il n'y a pas moyen d'employer la 
valeur 2 pour estimer le premier et le dernier trapèze, et la 
valeur f pour les autres, parce que l'on introduirait ainsi dans 
la formule les dérivées des ordonnées extrêmes. Poncelet a choisi 
une valeur intermédiaire, savoir 1, en posant 

1 
^ = ~(M ■+- m). 
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M. Parmentier, remarquant que le nombre des trapèzes inté- 
rieurs est en général supérieur à celui des extérieurs, a fait r = f , 

ce qui lui a donné 

1 

M •+- - m 

2 4 

p' -== —-_ (2M 4- m). 

2 

II. Nous sommes arrivésà la valeur tt^, plus rationnelle encore, 
de la manière suivante. Appelons u, la petite aire analogue à 
ANB comprise entre un arc ANB de la courbe, et sa corde AB 
dans un intervalle h (Fig. 9); l'aire ANBA 2 , comprise entre la 
tangente à une extrémité de Tare et la courbe, sera sensiblement 
2w, puisque la limite de leur rapport est égal à 2. L'aire analogue 
à ANB6 1 entre la courbe ANB et la demi-corde Ab { correspon- 
dant à un intervalle 2A, sera, pour la même raison, à peu près 
égale à bu. On aura donc approximativement 

M = S -t- 2n .2t* = S -4- inu, 
m = S — u — (2w — 2) hu — u = S — (8n — 6)u , 

On tire de là, en éliminant u, les relations approximatives 
(4n — 3) M -*- 2wm- = (6n — 3) S, 

S = v — - ==-(2M -h w) (M — 01)=*^, 

6n — 3 3 V ; 2» — 4 3 V ' 

en posant, comme au n° 7, 

1 14 4 4 4 

Wi = - (2M -«- m) --(M — m)=M M Ad. 

3 ' 2rc — 4 3 v ; 3 2rc — 4 3 

Les valeurs 2m, 4w, données plus haut, représentent exacte- 
ment les aires analogues à ANB, ANB6 1? dans le cas d'une para- 
bole du second degré. Il en résulte que la formule S =7r 1 donne 
exactement l'aire d'une parabole du second degré, et, par suite, 
d'après un principe général que nous établirons dans le para- 
graphe XVII, aussi celle d'une parabole du troisième degré. 



- 48 — 

III. Des considérations analogues à celles qui précèdent prou- 
vent que les valeurs/)",/)'" du n° 15 peuvent être remplacées 
par la valeur plus rationnelle . 

;r 2 ==l(2M'-*-m')-- — §(M ; — ro') 

= à (E + 21) 4- - faT - -hd". 

v ; 5 .2»-*-" 4 5 

La formule nouvelle, S = 7r 2 , est exacte aussi dans le cas d'une 
parabole du deuxième ou du troisième degré. 



1L\I 

Calcul de l'erreur par le théorème de Taylor. Formules 
de Simpson et de Weddle. 

36. Hypothèses préliminaires. Dans ce qui suit, nous allons 
ajouter les hypothèses suivantes à celles que nous avons faites 
jusqu'à présent : 1° Les dérivées F iv x, F v x, etc., qui vont entrer 
dans nos calculs ultérieurs, sont continues pour les valeurs de x 
considérées. 2° La dérivée d'ordre le plus élevé parmi celles que 
nous emploierons varie peu pour les valeurs de # considérées. 

L'erreur commise en employant les diverses formules données 
plus haut est de la forme 

e = aPOn) -*- 6F p (x 2 ) h- cF p (xz) 4- etc., 

F p x étant la dérivée d'ordre le plus élevé dont nous ayons besoin, 
a&j, # 2 , x$, ... des valeurs de x inconnues, mais comprises, entre 
les valeurs extrêmes de x. Cela posé, soit x^ une valeur de .t, 
comprise aussi entre ces valeurs extrêmes de x, on aura 

£ = (a + & + c + '") F p (x // ) •+- 
a[F*(à,) - F^(^)] 4- b[F>{x % ) - ¥ p {x F )] + c[F*{x z )— F*(a^)] + - 
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Par hypothèse, la partie de s écrite sur la seconde ligne est 
très petite. On pourra donc écrire approximativement 

f = (o + () + c + '")F p (x / ,). 

Cela revient à mettre, dès le début, dans l'expression de s, ou 
même dans les quantités qui ont servi à trouver s, x^ au lieu de 
ac 1? # 2 , x 3 , etc. C'est ce que nous ferons désormais, de sorte que 
toute égalité où entre x n'est qu'approximative. 

37. Formules de Simpson et de Weddle. I. L'élément de la 
première formule de Simpson, c'est-à-dire f h (y^ -h % 2 -+- y 2 ), 
devient, en posant j/. 2 = F (se), 

-h[i¥'(x) h- F'(x + h) -h F'(ar— ft)], 

et, en développant, 

2hF'+-h*F'" + — F'(ar„). 
3 36 v ^ y 



L'aire exacte est 



F(x + A) — F (a? — /i) = 2ÂF' -*- - /* 3 F'" + — F v (^.). 



\ 

3~~ 60 



L'erreur est donc approximativement 



36 l N 60 v ^ 90 . N 



Si F' (oc) est du second ou du troisième degré, F v (x)est nulle. 
Donc la formule de Simpson est exacte pour une parabole du 
deuxième ou du troisième degré. 

Si la base de l'aire à mesurer est 2nA, l'erreur pour chaque 
élément de base 2A est de la forme ~ W> F v (x^) et ces erreurs 
s'ajoutent, puisque F y (x) variant peu, ne peut pas, en général, 
changer de signe. 

4. 
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II. L'élément de la seconde formule de Simpson peut 
s'écrire 

-h\F'(x) +■ 3F'(x ■+• h) 4- 3F'(a? •+• 2A) -*- F'(x -+- 3A)] 
8 L 

9 9 27 33 

= 3AF' -+- - A 2 F"-t- - A 3 F'" -*- — A 4 F" + — /i 5 F v (x A 
2 2 8 16 v N 

L'aire exacte est 

9 9 97 81 

F(x -*-3A) — Fx=3AF'4- - A 2 F"-t- - A 3 F'"4- — A*F iv h- — -A 8 F v (x M ). 
2 2 8 40 



L'erreur est donc, approximativement, 



33 81 3 

; — W (x „) A 5 F V (ag = — i 



Cette erreur est nulle dans le cas d'une parabole du deuxième 
ou du troisième degré. 

III. L'élément de la formule de Weddle est : 

— A[F'(a; — 3A)-+- 5F'(x — 2A) + F'(x — A)-h6F'(x) 

-*- F'(x -4- A) -+- 5F' (oc + 2A) 4- F'(x 4- 3A)] 

81 7 

= 6F' -*- 9A 3 F'" 4- — A 8 F V 4- - A 7 F v "(x M ). 

L'aire exacte est 

81 243 

F(x -t- 3A) — F (x — 3A) = 6AF' 4- 9A 5 F'" 4- — F v 4- — /i 7 F v "(x^) . 

20 280 



Donc, 



7 243 1 

f = - A 7 F V " (x„) /a 7 F™ (xj = — A 7 F™(x A 

8 V ^ 280 V ^ 140 v N 



valeur nulle pour les paraboles du second, du troisième, du qua- 
trième et du cinquième degré. 
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IV. La valeur donnée par la seconde formule de Simpson 
a été établie de manière qu'elle fût exacte pour une parabole 
du troisième degré par Newton lui-même (Gauss, Werke, III, 
p. 202). 

La méthode qui a servi à trouver la formule de Weddle revient 
aussi à faire en sorte qu'elle donne l'aire exacte d'une parabole 
du 5 e degré dont on connaît sept ordonnées (voir plus bas, les 
notions historiques sur les diverses formules). Le second et le 
troisième résultat donnés plus haut, n'ont donc rien que de très 
naturel. 

Il n'en est pas de même du premier. La valeur approchée 
1 h(jj) 4- 4t/ 2 -+- Vz) de l'élément a été obtenue, comme nous le 
dirons dans l'historique, en cherchant l'expression qui donne 
l'aire exacte d'une parabole du second degré. Cependant, cette 
formule est exacte aussi, comme nous venons de le voir, pour 
une parabole du troisième degré. C'est là un théorème curieux (*) 
dont la formule de Gauss exprime, au fond, une belle générali- 
sation pour les paraboles d'ordre supérieur à ordonnées non 
équidistantes (Gauss, Werke, III, pp. 165-196,202-206). 

Il en existe une autre généralisation pour les paraboles quel- 
conques d'ordre pair, comme nous allons le montrer, avant de 
continuer la recherche de l'erreur pour les autres formules. 



Théorèmes sur les paraboles. 

38. Théorème I. Soit 

h (ay ï h- by* -*- cy z -t- . . . -*- gy n ), 
une formule donnant Vaire exacte d'une parabole d'ordre p, p 



O Ce théorème, donné par M. Catalan, dans les Nouvelles Annales de Mathématiques, 
i Te série, t. XVI, p. 312, et dans le Manuel des Candidats à l'école polytechnique, t. II, p. 
295, a été retrouvé par M. Parmentier {Congrès de Montpellier de l'Association française 
pour l'avancement des Sciences, 1879). 
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étant inférieur à ri, au moyen de n ordonnées de cette parabole , 
distantes entre elles d'une quantité h. Si on applique cette formule 
à une courbe y = F'x, V erreur commise s sera de la forme 

En effet, considérons une parabole d'ordre p, y = f (oc), 
osculatrice à y = F' (oc), en un point quelconque [X, F'(X)], 
de sorte que 

f(X)-F'(X), /"(X) = F"(X), ,.., p^(X) = F^(X). (1) 
Soient 

x 1 = X -4- M, x 2 = # 1 •+- />, x 3 = a; 2 H-/i ? ... tf n = x n _ 4 -t- A , 
lès abcisses correspondant à y { , i/ 2 , 2/ 3 , ..., ?/„. Posons : 
S = F(ag-F(* 4 ), 

s = fM — f( x i)> 

S' = h [a¥ f (x { ) -+- 6F' (x 2 ) -+- • • • -t- flfF' (*?„)], 
«' = h [af (x â ) + 6/' (x a ) -*- . . . -h gf (x„)], 

Développons, suivant les puissances de A, la différence S — s, 
en ne laissant dans le développement d'autre abcisse que X. A 
cause des relations (1), on aura évidemment 

S — s = A i hr**F"*(x /l ) i 

car les autres dérivées de F, f sont égales et / p+ " 2 (x) = 0. Pour 

la même raison, 

S' — s r = A 2 ^+ 2 F^ 2 (^), 

Mais, d'après Phypothèse faite sur la formule qui donne s', on a 
s' == s. Donc : 

£ = S — S' = (A* — A 2 ) A" +a F" +2 (^) = AA^F^ 2 ^); 

ce qu'il fallait démontrer. 
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En particulier, si p = (n — 1), on a S — S' = A h n+1 F (a^), 
ce qui est la justification de la méthode de quadrature approchée, 
due à Newton, dont il est parlé au n° 1, pour les courbes dont 
1 ordonnée F' (x) vérifie les conditions énoncées au n° 36 (*). 

39. Théorème II. Soit 

h [a [y, n- y n ) -*- 6 (y 2 4- y n ^) + c{y z + y n .^) -+■•..], 

tme expression donnant Faire exacte d'une parabole d'ordre 2q , 
2q éftm/ inférieur à n, au moyen de n ordonnées de cette para- 
bole, cette expression étant symétrique par rapport aux ordonnées 
également éloignées des extrêmes. Si on applique cette formule à 
une parabole d'ordre 2q -4- \, elle en donnera aussi Faire exacte. 
Soient, en effet, y = F' (x), l'équation de cette parabole, et 
2X = Xi H- x n . On aura 

valeur qui, développée , ne contient que des puissances impaires 
de A. Ensuite : 

a 

i r / i . ' ■_. ,\ _. / i . . m 

etc.] 



F'(x-i(fi-l)A)H-F'(x + i(n-*)A)J. 
[6 F' (x — -(n — 5)/*) + F' (x + * (» — 5) a)1 ■ 



valeur qui, développée, ne contient aussi que des puissances im- 
paires de h. D'après le n° précédent, le premier terme de la dif- 
férence S = S' est 

ou plutôt, 

AA^F'^X),. 



(*) Si l'on veut justifier la méthode de Newton, en employant des tracés graphiques, on 
reconnaît bien vite qu'à partir du troisième degré, très souvent les paraboles s'écartent 
beaucoup des courbes avec lesquelles elles ont autant de points communs que possible, 
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puisque F' (oc) est d ordre (2g -+- 1) et F 234 " x constant. D'autre 
part, cette différence ne contient que des puissances impaires 
de h; donc A = 0, c'est-à-dire que S = S'. 
Corollaire. Si l'on applique la formule 

h [a (y t -h y n ) + 6 (y, h- y n _ { ) + •••], 

d wne courbe y = F' (x), l'erreur sera de la forme 

A/^+*fw ( Xfl ). 

XVIII 

Calcul de l'erreur. Formule des trapèzes, de Poncelet, etc. 

40. Formules des trapèzes, de Poncelet, de M. Parmentier, 
de M. Dupain et formule nouvelle. I. D'après ce que Ton a vu au 
n° 32, Terreur que comporte la formule des trapèzes est égale, 
approximativement, à 

(I-i)^"(x,)=l^"(^). 

IL Dans la formule de Poncelet, Terreur, pour les éléments 
extrêmes, est 

*_i)*V»^-l*«P'»(x,). 
Pour les éléments intérieurs, on a, d'après le même numéro, 

III. Dans la formule de M. Parmentier, Terreur, pour les élé- 
ments extrêmes, est 

i\ I \ 1 
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Pour les éléments extérieurs, Terreur, au contraire, s'exprime 
en h* 1 , comme il est évident, d'après 32, II et 36. Si Ton prenait 
l'erreur totale pour tous les éléments, on pourrait faire en sorte 
que Terreur ne contienne que des puissances impaires de /*, 
comme on Ta vu au n° 39. L'erreur dans la formule de M. Par- 
mentier est donc du troisième degré en h, mais cette erreur en A 5 
provient uniquement des éléments extrêmes. 

IV. Enfin les formules S ===== 7^, S = tt 2 étant exactes pour des 
paraboles du deuxième, et, par suite (n°39),du troisième degré, 
l'expression de Terreur est de la forme 

£ = A/fc 5 F v (^). 

41. Formule de M. Catalan. L élément cr dont il est parlé 
au n° 26, a pour aire exacte 

F(x + h) — ¥x = hY' + -h 2 ¥" H-ifc'F'"'-*-— * 4 P f (a?A 

\ / . ^ 6 24 p' . 

La valeur approximative est : 

\h[W{x) + W(x + h) — F'(a;-t-2À)] 

1 À 



4 1 

2 6 



: hF •+- z Ii 2 F" h- - A 3 F'" -+- O.A'F"^). 



L'erreur est donc 

1 

24 

nulle dans le cas d'une parabole du second degré, et, par suite 
du troisième degré. On peut déduire ce résultat aussi du n° 32, 
III, et de la théorie de la parabole. 

Si Ton prend pour F' (oc) l'ordonnée moyenne de Taire consi- 
dérée, il est clair, d'après le n° 39, que Terreur sera de la forme 

AtfF'(^), 

comme pour les deux formules nouvelles et pour celles de 
Simpson. 
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Comparaison des formules. 

42. Comparaison an point de vue de V exactitude. Si l'on géné- 
ralise les résultats trouvés dans le § XII, en traitant divers exem- 
ples, et ceux que nous avons obtenus dans les §§ XVI et' XVIII, 
moyennant les hypothèses du n° 56, on peut diviser les for- 
mules étudiées en trois groupes. 

Le premier comprend les deux formules de Simpson et celle 
de Weddle, qui sont exactes, les deux premières, pour les para- 
boles du 2 me ou du 5 me degré, la dernière pour celles du 2 me , 
du 3 me , du 4 me et du 5 me degré. L'erreur pour les deux pre- 
mières est de la forme AhW y (x fl ) 9 pour la dernière de la forme 
BA 7 F V,I (^). Pour arriver à cette estimation de Terreur, on sup- 
pose, pour les premières, que F'(x), F"0*0> F'"(ac), F IV (V)> F v (x) 
sont continues et que F y (x) varie peu, dans un intervalle 2A pour 
la première, dans un intervalle 5/& pour la seconde. Pour la for- 
mule de Weddle, on admet que F'(x) et ses six premières déri- 
vées sont continues dans un intervalle 6h et que F v "(œ) varie 
peu dans le même intervalle (*). Â cause de ces hypothèses plus 
étendues et plus nombreuses, il est probable que la formule de 
Weddle ne peut pas être mise au-dessus de celles de Simpson. 
D'ailleurs, quand on compare les formules ici en question, d'après 
les limites supérieures de Terreur commise, limites trouvées 
§§ IV, X, XI, on reconnaît que la première formule de Simpson 
l'emporte sur les autres. On peut donc, en réunissant les divers 
éléments d'appréciation, conclure que les trois formules du pre- 
mier groupe ont à peu près la même valeur, au point de vue de 
l'exactitude. 



(*) Une formule obtenue par la méthode de Côtes et exacte pour une parabole d'ordre n, 
conduit à une erreur de la forme. A/i n + 2 F«-t- 2 (#/*), pour une courbe t/ = F(x), moyennant 
(rc-hl) conditions de continuité, dans un intervalle nh, et pourvu que F n + i {x) varie peu 
dans cet intervalle. Le nombre des conditions pour que la formule soit très exacte croît 
donc très rapidement avec le degré d'exactitude. 
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Le second groupe contient les formules de M. Catalan et de 
M. Dupain et la formule nouvelle. Elles sont exactes aussi pour 
les paraboles du 2 me et du 3 me degré; l'erreur est de la forme 
AhPFXxp); mais pour arriver à ce résultat, nous avons supposé 
cfue F'(x) et ses quatre premières dérivées sont continues dans 
un intervalle nh, où d'ailleurs F v (x) varie peu. A cause de la 
grandeur de cet intervalle, en général, ces trois formules sont 
probablement moins exactes que les précédentes. 

Le troisième groupe comprend les autres formules. Celle de 
M. Parmentier donne un résultat très voisin de celui que fournit 
la formule de M. Dupain; l'erreur ne dépend de A 3 qu'à cause 
des éléments extrêmes. Il est donc extrêmement probable qu'elle 
vaut mieux, en général, que celle de Poncelet, et à fortiori que 
la formule des trapèzes (*).. 

43. Comparaison au point de vue de la simplicité des calculs. 
Au point de vue de la simplicité des calculs, les neuf formules 
considérées forment encore trois groupes, presque identiques aux 
précédents. 

Le premier comprend la formule des trapèzes, qui est la plus 
facile à retenir, puis les formules de Poncelet et de M. Parmen- 
tier qui contiennent moins d'ordonnées, enfin la formule de 
M. Dupain, qui est la moins simple des quatre. 

La formule de M. Catalan et la formule nouvelle forment le 
second groupe : elles sont un peu plus compliquées, à cause de 
la quantité d! qui y entre. 

Enfin, les formules dont l'application entraîne les calculs les 
plus pénibles sont celles de Simpson et de Weddle, qui consti- 
tuent le troisième groupe. 

44. Conclusion. Si l'on n'a pas besoin d'une grande approxi- 
mation, le mieux, d'après ce qui précède, est de recourir à la 



On trouve dans la Statique graphique de Culmann, traduction française, Paris, 
Dunod,1880,pp. 403-105, une discussion sur la valeur comparée des formules de Simpson 
et de M. Parmentier où il essaie de prouver que la dernière est supérieure à la première. 
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formule des trapèzes, ou, si le nombre des ordonnées est impair, 
à celle de M. Parmentier, qui est plus exacte et exige le calcul 
de beaucoup moins d'ordonnées. 

Si Ton veut une exactitude plus grande, et si le nombre des 
ordonnées est pair, on prendra la formule de M. Catalan. Quand 
le nombre des ordonnées est impair, on se servira plutôt de la 
formule de M. Dupain, qui emploie moins d'ordonnées, ou, 
si l'exactitude est requise avant tout, de la première formule de 
Simpson. 

Les autres formules peuvent être employées dans des cas spé- 
ciaux, par exemple, la seconde formule de Simpson, dans le cas 
où la base de l'aire à mesurer est divisée en 3n parties égales. 



Historique. 

45. Formules de Simpson et de Weddle. Thomas Simpson (*) 
a publié sa formule dans ses Mathematical Dissertations, 1743, 
p. 109 (**). Au fond, il cherche l'aire des trapèzes curvilignes 
analogues à aABCc(Fig. 10), en substituante la courbe ABC une 
parabole passant par les trois points A, B, C et dont l'axe est per- 
pendiculaire à ac. Cette démonstration a passé dans la plupart 
des traités de calcul intégral. — Poncelet, dans son Introduction 
à la Mécanique industrielle (n° 180, pp. 187-190 de la 2 me édi- 
tion, Metz, 1839; 3 me édition, Paris, 1870, pp. 194-197), en a 
donné une démonstration élémentaire comme il suit : Pour trou- 



(*) Né en 4710, mort en 1761; ne pas le confondre avec son contemporain, Robert 
Simson, né en 1687, mort en 1761 aussi. 

[**) Nous empruntons ce renseignement à une note de la page 67, de la célèbre édition 
de la Résistance des Corps solides de Navier, donnée par M. B. de Saint -Venant, Paris, 
Dunod, 1864. T. Simpson, dans son ouvrage intitulé : The Doctrine and Application of 
Fluxions, London 1750, p. 202, fait allusion à sa formule, mais ne la donne pas. Probable- 
ment, les Dissertations de Simpson contiennent aussi sa seconde formule, que nous avons 
empruntée à un ouvrage de Rankine [Useful Rides and Tables. Fourth Edition. London, 
C.GriffinandO, 4873, p. 68). 
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ver Taire aABCc, il divise ac en trois parties égales aw, uv 9 vc, 
et mène les ordonnées de la courbe wU, vY, ainsi que les cordes 
AU, UV, VC, dont la seconde UV rencontre 6B en W. L'aire 
curviligne aABCc est un peu plus grande que celle du polygone 
aAUVCc qui a pour mesure ~ h (y^ -h 46W > -+- j/ 3 ) ; en rempla- 
çant 46 W par la quantité un peu plus grande ây 2 , ce qui revient 
à ajouter au polygone aAUVCc, quatre fois le petit triangle UBV, 
on se rapprochera vraisemblablement de l'aire exacte de aABCc. 
— MM. Rouché et de Comberousse, dans leur Traité de Géomé- 
trie (I re partie, Paris, Gauthier-Villars, 1879, pp. 294-296) ont 
simplifié cette démonstration en remarquant que si AC coupe 6B 
en 6 1? Faire du trapèze a A Ce est égale à | h (^ -t- 466 d -+■ y z ) ; 
remplaçant kbb i par 4ï/ 2 > ce qui revient à ajouter |au trapèze les 
| du triangle ABC, on se rapprochera, vraisemblablement encore, 
de Taire cherchée. — Il est clair qu'aucune de ces trois méthodes 
ne fournit une vraie démonstration de la formule, puisque Ter- 
reur maxima commise n'est pas évaluée (*). 




Fig. 10. 



(*) Dans les Nouvelles Annales de Mathématiques, l re série, t. XV, pp. 291-293, un 
Anonyme metjla formule de Simpson sous la forme s = t 4- 1 {t — m') et dit que Saygey 
en a donné une démonstration simple dans sa Géométrie élémentaire. Nous n'avons pu 
consulter cet ouvrage. 
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Nous ignorons où la formule de Weddle a été publiée pour la 
première fois. Nous Pavons empruntée à la traduction allemande 
du Calcul des différences de Boole, qui l'obtient en assimilant la 
courbe dont on cherche Taire à une parabole du 6 me degré et 
altérant légèrement le résultat, de manière qu'il soit exact pour 
une parabole du 5 me degré (*). 

4©. Formule de Poncelet, de M. Parmentier et de M. Dupain. 
La formule de Poncelet, exposée par lui dans son Cours de la 
faculté des sciences de Paris, n'a pas été publiée par son auteur, 
mais elle s'est rapidement répandue. D'après M. de Saint-Venant, 
dès 1851, elle est donnée déjà dans le Traité élémentaire de Méca- 
nique de ML Sonnet. On la trouve maintenant dans la troisième 
édition de la Mécanique industrielle, publiée par M. Kretz, 
pp. 197-201. 

La formule de M. Parmentier a été trouvée, en 1854, par ce 
savant officier (alors capitaine, maintenant général), et publiée, 
dans le n° 16 du Mémorial de V officier du Génie vers le milieu 
de cette année, puis dans les Nouvelles Annales de Mathématiques, 
l re série, t. XIV (1855), pp. 570-384. La démonstration, au fond, 
ne diffère pas de celle qui a été donnée plus haut. Le général 
Piobert arrivait, quelques mois après, à la même formule, par une 
voie moins nette, dans un article publié en septembre 1854, dans 
les Nouvelles Annales de Mathématiques, l re série, t. XIII, pp. 525- 
551. M. le général Parmentier est revenu sur sa formule, en 1876, 
dans un article intéressant, inséré dans le tome XV, 2 e série 
pp. 241-251, du même recueil, où il donne, comme nous l'avons 
fait plus haut, l'estimation de l'erreur maxima qu'elle comporte. 
Il remarque, de plus, qu'on arrive à cette formule en retranchant 
de 2s (Fig. 2) l'aire des trapèzes aABb, kKLl, puis l'aire curvi- 
ligne bBIikb, calculée par la formule de Simpson. Cette nouvelle 
démonstration de la formule de Parmentier fait saisir, sous un 
nouveau point de vue, la raison de son exactitude supérieure à 



(*) Die Grundlehren der endlichen Differenzen- und Summenrechnung von 
C. Boole. Deutsch bearbeitet von C.-H. Schnuse. Braunschweig, Leibrok, 4867, pp. 40-41. 
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celle de la formule de Poncelet (*). Un travail un peu plus étendu 
sur le même sujet, avait été présenté parM.Parmentierà l'Associa- 
tion française pour l'avancement des sciences, dans sa réunion de 
Nantes, en 1 875, et se trouve dans le Compte rendu de cette réunion, 
pp. 179-200. Ce sont ces divers mémoires de M. Parmentier qui 
nous ont inspiré les recherches contenues dans le présent travail. 
La formule de M. Ch. Dupain, qui a passé inaperçue, malgré 
sa valeur théorique et pratique, a été publiée en 1 858, dans les 
Nouvelles Annales de Mathématiques (t. XVII, pp. 288-295). 
L'auteur l'obtient en prenant une expression ttj de la forme 
^ ™ ? te " e 9 ue l' erreur S — 7T 4 soit d'ordre le plus élevé pos- 
sible en h. Il lui attribue la même exactitude qu'aux formules de 
M. Catalan et de Simpson. Il calcule à peu près comme nous, 
par le théorème de Taylor, la valeur de l'erreur que comportent 
les formules de Poncelet, de M. Parmentier, de M. Catalan et de 
Simpson, mais chose curieuse, il semble ne connaître la limite 
supérieure de l'erreur que pour celle de Poncelet. 

4*. Formule de M. Catalan. M. Catalan a publié sa formule 
en 1851 (Nouvelles Annales de Mathématiques, I re série, t. X, 
pp. 412-415). Il l'obtient en regardant l'aire curviligne aABô, 
comme celle d'une parabole à axe parallèle à Aa, et passant par 
A, B, C, l'aire curviligne 6BCc, comme celle d'une parabole ana- 
logue passant par B, C, D; et ainsi de suite jusque l'aire ilKLl 
qu'il calcule comme Simpson. Il refait ensuite des calculs ana- 
logues en sens inverse, en partant de L et marchant vers A, de 
sorte que l'aire ILKk est maintenant regardée comme celle d'une 
parabole analogue passant par L, K, I et ayant un axe parallèle 
à Aa; kKli, l'aire d'une parabole analogue passant par K, I, L; 
enfin l'aire cCBAa est calculée comme chez Simpson. La moyenne 
des deux résultats obtenus donne la formule de M. Catalan. 



(*) Cette exactitude supérieure a été reconnue par Poncelet lui-même. Malgré cela, 
M. Kretz, dans la troisième édition de la Mécanique industrielle, en se basant sur la note 
de Dupain analysée plus bas, dit que la formule de Poncelet est préférable à celle de Par- 
mentier, assertion vraiment regrettable, parce qu'elle se trouve dans un livre classique. 
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NOTE 



Une formule inédite de M. Parmentier. 

M. le général Parmentier nous autorise à publier ici une formule 
inédite, à la fois très simple et très exacte, qu'il a trouvée autrefois 
en cherchant à perfectionner celle que nous avons donnée au n° 7. 

Cette formule peut s'écrire, en se servant des notations de notre 

mémoire, 

S ==£>! (approximativement), 

1 1 1 

p. = 2/iP M h — hd' = p-*- - hd'. 

™ 2 6 ■ 6 

Il est clair que p l5 étant inférieur à 2AP ou M ((/surpasse d% 
d'après le n° 10), et supérieur à p, et, par suite, à m, Terreur que 
comporte la nouvelle formule est inférieure à la plus grande des deux 
différences M — pi,Pi — m, c'est-à-dire à la seconde, ou à ihd-+-$hd', 
et, à fortiori, à § hd. 

On reconnaît aisément que la formule S=p 4 est exacte pour une 
parabole du 2 e ou du 5 e degré, par des raisonnements analogues à 
ceux du n° 35 et par le théorème du n° 39. 

Si on l'applique aux quatre exemples du § XII, on trouve : 

Pi=S; p^S — 65; p, = S-i-89; p, = S-Ml. 

Enfin, on peut écrire 

p t î=3s — 2C, 

relation inattendue entre la nouvelle formule et celles de Simpson et 
de M. Catalan. On peut en conclure qu'il est inutile d'employer cette 
dernière dans le cas où la base de l'aire à mesurer est divisée en un 
nombre pair de parties égales, puisque la nouvelle formule, qui 
comporte des calculs beaucoup plus simples, est probablement tout 
aussi exacte (*). 

(*) Ce Mémoire est extrait des Annales de la Société scientifique de Bruxelles, t. V, 
pp. 231-290. Dans les Annales, la formule S = 7T t est donnée comme une formule nou- 
velle, parce que nous ignorions l'existence de la petite note de M. Dupain, qui nous a été 
signalée par M. Parmentier. Nous saisissons cette occasion pour remercier le savant 
général des précieux renseignements qu'il a bien voulu nous communiquer sur la ques- 
tion des quadratures approchées. 
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